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ABSTRACT 

We derive metal abundance estimates for a large sample of starbursting emission-line galaxies (ELGs). 
Our sample is drawn from the KPNO International Spectroscopic Survey (KISS) which has discov- 
ered over 2000 ELG candidates to date. Follow-up optical spectra have been obtained for ~ 900 of 
these objects. A three step process is used to obtain metal abundances for these galaxies. We first 
calculate accurate nebular abundances for 12 galaxies whose spectra cover the full optical region from 
[OII]AA3726,29 to beyond [SII]AA6717,31 and include detection of [OIH1A4363. Using secondary metal- 
licity indicators R23 and p^, we calculate metallicities for an additional 59 galaxies with spectra that 
cover a similar wavelength range but lack [OIIIJA4363. The results are used to calibrate relations between 
metallicity and two readily observed emission-line ratios, which allow us to estimate coarse metallicities 
for 519 galaxies in total. The uncertainty in these latter abundance estimates is 0.16 dex. From the 
large, homogeneously observed sample of star-forming galaxies we identify low metallicity candidates for 
future study and investigate the metallicity-luminosity relation. We find a linear metallicity-luminosity 
relation of the following form: 12+log(0/H)= 4.059 - 0.240M B , with an RMS scatter of 0.252. This 
result implies that the slope of the metallicity-luminosity relation is steeper than when dwarf galaxies 
are considered alone, and may be evidence that the relationship is not linear over the full luminosity 
range of the sample. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: starburst 



1. INTRODUCTION 

Observing metal abundances in galaxies is a valuable 
probe of galaxian star-formation histories and chemical 
evolution. Assuming that all galaxies begin with the same 
primordial abundances of elements, roughly 75% hydrogen 
and 25% helium, measurements of heavier elements indi- 
cate subsequent star formation accompanied by supernova 
explosions which enrich their surroundings with metals. 
Recent studies show that metallicity is linked with galaxy 
luminosity (Skillman et al. 1989; Zaritsky et al. 1994; 
Richer and McCall 1995; Garnett et al. 1997; Hunter and 
Hoffman 1999; Pilyugin and Ferrini 2000; Pilyugin 2001b) 
in that more luminous galaxies tend to be more metal rich 
than less luminous galaxies. This may indicate an evo- 
lutionary trend with several possible explanations. For 
example, a simple closed box model gives rise to such a 
trend (Hidalgo-Gamez and Olofsson 1998). Alternatively, 
evidence for large disruptions in the gas content of galax- 
ies from supernova explosions may play a significant role. 
Specifically, supernovae have been proposed as a mech- 
anism for removing large amounts of metal-enriched gas 
from low-mass systems (e.g., Mac Low and Ferrara 1999). 
This in turn may explain the relative dearth of metals 
in dwarf galaxies. In addition, higher astration levels in 
more luminous galaxies may contribute to their increased 
metallicity per mass (Pilyugin and Ferrini 2000). 

The KPNO International Spectroscopic Survey (KISS) 
has identified over 2000 emission-line galaxy (ELG) candi- 
dates ranging in absolute magnitude from Mb = -22 to -12 
(Salzer et al. 2000, 2001). As such, it provides a large sam- 



ple of galaxies for which metallicities can be derived and 
from which the metallicity-luminosity relation can be stud- 
ied. The survey lists include massive starburst-nucleus 
galaxies, intermediate-mass irregular galaxies, low-mass 
dwarf irregulars and blue compact dwarfs. In this work we 
derive coarse metallicity estimates for the 519 starburst 
galaxies which have adequate quality follow-up spectra. 
The results provide a list of low metallicity candidates for 
future study as well as the largest sample of galaxies to 
date for use in studying the metallicity-luminosity relation 
Most previous studies of the metallicity-luminosity re- 
lation have concentrated on low metallicity galaxies where 
accurate abundances are readily available. Studies of ir- 
regulars with current star formation (Skillman et al. 1989; 
Richer and McCall 1995) show evidence for a linear re- 
lation. However, another recent study of dwarf irregu- 
lars (Hidalgo-Gamez and Olofsson 1998) does not support 
these results. In re-examining the Hidalgo-Gamez and 
Olofsson data, Pilyugin (2001b) believes that noise in the 
[OIIIJA4363 line is responsible for the lack of a relation in 
their sample. 

On the high mass end, Garnett et al. (1997) compiled a 
data set of 29 luminous spirals, with metallicities derived 
from HII regions. The data set is taken from Vila-Costas 
and Edmunds (1992), Zaritsky et al. (1994) and Ryder 
(1995). They find that these more massive galaxies follow 
a similar trend to dwarf irregulars of increasing metallicity 
with luminosity. Garnett et al. do not offer a specific fit to 
the relation, but rather state that the high luminosity re- 
sult maps smoothly onto the dwarf irregular relations such 
as that of Skillman et al. The same conclusions are drawn 
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by Zaritsky et al. (1994). Pilyugin and Ferrini (2000) 
combine the low and high mass galaxy samples to provide 
a new fit to the metallicity-luminosity relation and show 
a somewhat steeper relation than Skillrnan ct al (1989). 

By using the large KISS sample of galaxies we investi- 
gate whether the metallicity-luminosity relation is actually 
present on all mass, luminosity and mctallicity scales. Our 
abundance determinations are coarse and we use all mor- 
phological types for our metallicity-luminosity relation. 
Therefore our relation contains significantly more scatter 
than previous results. However, our luminosity range goes 
three magnitudes brighter than most previous studies, and 
we use up to 20 times more galaxies, giving us a better view 
of the overall form of the metallicity-luminosity relation, 
as well as its intrinsic scatter. We confirm the existence 
of a metallicity-luminosity relation that varies smoothly 
from massive luminous galaxies, Mb < —21, to blue com- 
pact dwarfs, Ms > —16. However, when the metallicity- 
luminosity relationship is calculated using the wide range 
of galaxy types found in the KISS sample, the slope of 
the relation is steeper than indicated by the studies of 
dwarf galaxies such as Skillrnan et al. and Richer and Mc- 
Call. As a result we find that a simple extrapolation of 
the dwarf galaxy relationship to more massive systems is 
inappropriate. This may be an indication that the overall 
form of the metallicity-luminosity relation is not linear but 
rather requires a higher order polynomial to fit the data. 
It may also indicate that different galaxy types obey dif- 
ferent metallicity-luminosity relations. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations 

Drawing from the KISS sample of emission-line galax- 
ies, we gathered both imaging and spectral data for ^900 
galaxies over a four year period from 1998 - 2001. Pho- 
tometry is presented in the survey lists (Salzer et al. 2001; 
2002a; Gronwall et al. 2002a) while results of spectral 
follow-up are given in a series of papers (Melbourne et al. 
2002, hereafter Paper II; Salzer et al. 2002b; Wegner et 
al. 2002). A summary of the overall spectroscopic prop- 
erties of the KISS ELGs will be presented in Gronwall et 
al. (2002b). 

The spectral data can be classified into 3 groups based 
on spectral coverage and quality. Group I spectra, an ex- 
ample of which is given in Figure la, cover the full opti- 
cal region from [OII]AA3726,29 to beyond [SII]AA6717,31 
and contain the [OIIIJA4363 line necessary for accurate 
abundance measurements. Abundances are calculated for 
Group I galaxies using the [OIIIJA4363 line as an electron 
temperature indicator. All 12 objects in this group were 
observed with the Lick 3m telescope. The spectroscopic 
data and abundance analysis are presented in detail in Pa- 
per II. 

Group II data also cover the full optical region from 
[OII]AA3726,29 to beyond [SII]AA6717,31, but do not nec- 
essarily contain the [OIII]A4363 line needed for accurate 
abundance measurements. Metal abundance estimates are 
derived from the strong oxygen lines [OII]AA3726,29 and 
[OIII]AA4959,5007 using the secondary metallicity indica- 
tors i? 23 (Pagel et al. 1979) and p 3 (Pilyugin 2000). The 
data were obtained at the Lick 3m, the APO 3.5m and 
KPNO 2.1m telescopes and include 59 additional objects. 



A sample spectrum is shown in Figure lb. 

Group III spectra do not reach blueward to the 
[OII]AA3726,29 line. Metal abundances are derived 
for these galaxies based on empirical relations between 
the [OIIIJA5007/H/3 line ratio and mctallicity and the 
[NII]A6583/Ha line ratio and metallicity. Data in this 
category were obtained from the above mentioned tele- 
scopes as well as the WIYN 3.5m, the MDM 2.4m and 
the Hobbey-Eberly 9m (HET). As the HET is a 9m-class 
telescope, the spectra have a high signal-to-noise ratio as 
shown in Figure lc. Unfortunately, they do not extend 
much bluer than Hy. The [SII] doublet on the red side 
is also often beyond the wavelength range of the spectro- 
graph. The MDM images have a low dispersion. There- 
fore, as is shown in Figure Id, the [NIIJA6583 line blends 
with Ha, and the [SII]AA6717,31 doublet is also blended. 

Details of the observations, data reduction methods and 
line measurements are given in a series of papers which 
present the results of our extensive spectroscopic follow- 
up (Paper II; Salzer et al. 2002b; Wegner et al. 2002). 

3. METALLICITIES OF KISS GALAXIES 

The most accurate calculation of metallicity in nebu- 
lar star- forming regions uses the [OIIIJA4363 line to mea- 
sure the electron gas temperature (Osterbrock 1989; Izo- 
tov et al. 1994). This method, referred to here as the 
T e method, works well for low metallicity systems (gener- 
ally 12 + log(0/H) < 8.2) where the [OIII]A4363 line is 
observable. For systems of higher metallicity, [OIII]A4363 
is often too weak to observe or too noisy to be trusted. 
However, the strong nebular lines alone contain the nec- 
essary information to arrive at relatively good estimates 
of the oxygen abundances in star-forming regions (Pagel 
et al. 1979; McGaugh 1991; Pilyugin 2000). The tradi- 
tional i?23 method (Edmunds and Pagel 1984) typically 
results in metallicities within 0.2 dex of the more exact 
T e abundances. Pilyugin (2000) improves on the standard 
i?23 method by introducing his p^ factor calculated from 
the strong oxygen lines. The p^ factor replaces the tem- 
perature as a descriptor of the conditions in the nebula 
and allows his method to correlate with the T e method to 
within 0.1 dex for starbursts with metallicities below 12 + 
log(0/H) < 7.9. In this section we make use of the T e , 
i?23 and P3 methods to estimate metallicities for 71 KISS 
galaxies. These data are then used to show that both 
the [OIII]A5007/H/3 and [NII]A6583/Ha line ratios corre- 
late with metal abundance. We use the relationships to 
estimate coarse metallicities for the large, homogeneously 
observed KISS sample of galaxies. The details of each step 
follow. 

3.1. T e Metallicities 

Spectra of 12 KISS galaxies contain the necessary infor- 
mation for high quality metal abundance determinations. 
Following the standard procedure (Osterbrock 1989; Izo- 
tov ct al. 1994) we calculate the electron density from 
the [SII] line ratio, and the electron temperature from the 
[OIII] line ratio. Metal abundances are derived using the 
IRAF NEBULAR package (de Robertis et al. 1987; Shaw 
et al. 1995). Details of the observations, data reductions 
and analysis are given in Paper II. The results are pre- 
sented in Table 1. We use the T e abundance results in 
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Fig. 1. — Sample spectra from the KISS archive: a) A Group I spectrum taken at the Lick 3m telescope. It covers the full optical region 
from [Oil] AA3726,29 to beyond [SII] AA6717,31 and includes detection of the [OIII]A4363 line necessary for accurate metal abundances, b) A 
Group II spectrum taken at the Kitt Peak 2.1m telescope. It also covers the full optical region but does not have the signal-to-noise ratio 
necessary for the detection of [OIIIJA4363. c) A Group III spectrum taken with the HET 9m telescope. It has a high signal-to-noise ratio 
because of the large telescope aperture, but does not cover the full optical region and is specifically lacking [Oil] AA3726,29. d) Also a group 
III spectrum, taken with the MDM 2.4m. It has lower dispersion; the nitrogen lines are blended with Ha, and the sulfur doublet is also 
blended. 
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Fig. 2. — This line diagnostic diagram shows the relationship between R23 and the [Nil] A6583/Ha line ratio. Metallicity varies smoothly 
along the observed distribution, with high metallicity galaxies to the lower right and low metallicity galaxies to the upper left. When 
log([NII]A6583/Ha) < -1.3 the galaxy is on the low metallicity branch of the R23 parameter. When - 1.3 < log([NII] A6583/Hct) < -1.0 the 
galaxy is in the turn around region, and when log([NII] A6583/Ha) > -1.0 the galaxy is the high metallicity branch of the R23 parameter. 



Section 3.3 when we correlate metallicity with emission- 
line ratios. 

3.2. Secondary Metallicity Indicators R23 and p 3 

The strong oxygen lines [OIII]AA4959,5007 and 
[OII]AA3726,29 contain the necessary information to pre- 
dict the metallicity of an HII region (Pagcl ct al. 1979; 
McGaugh 1991). Traditionally this has been done by in- 
voking the i?23 parameter (Pagel et al. 1979), where; 

/([OIII]A4959 + A5007) + /([OII]A3726 + A3729) 



R 



23 



/(H/3) 



(1) 

The i?23 parameter has been correlated with metallicity by 
measuring oxygen abundances via the T e method for large 
samples of galaxies and HII regions. A complication with 
this method is that the dependence of metallicity on R23 is 
double valued. On the low metallicity end, 12 + log(0/H) 
< 7.9, i?23 increases with metallicity. As the amount of 
oxygen in the nebula increases the strength of the emis- 
sion leaving the nebula in the forbidden oxygen lines also 
increases. For higher metallicities, 12 + log(0/H) > 8.1, 
i?23 decreases with increasing metallicity. In this case the 
majority of the energy leaves the nebula in other nebular 
lines. Between the high and low metallicity branches is the 
"turn-around region" where R23 ceases to be a good pre- 
dictor of metallicity. In the "turn-around region" galaxies 
with the same -R23 ratio can have a fairly wide range of 
metallicities. 

We distinguish between the upper and lower branches 
of the i?23 relation by observing the [NII]A6583/Ha line 
ratio. This is demonstrated in Figure 2, a line diagnos- 



tic diagram plotting log(i?23) vs. log([NII]A6583/Ha) 
. In this plot, metallicity varies smoothly across the fig- 
ure with low abundance galaxies found in the upper left 
hand corner and high abundance galaxies in the lower 
right. We adopt the following criteria for determining 
on which branch of the i?23 relation a given spectrum 
lies. We assign objects with log([NII]A6583/Ha) < -1.3 
to the low metallicity branch of the i?23 relation. Ob- 
jects with log([NII]A6583/Ha) > —1.0 are associated 
with the high metallicity branch. Objects with -1.3 < 
log([NII]A6583/Ha) < -1.0 are considered turn around re- 
gion objects. We do not calculate i?23 metallicities for 
these galaxies. 

Unfortunately R23 does not correlate perfectly with 
metallicity Galaxies with a given value for i?23 can exhibit 
a range of T e -method-determined metallicities. Other fac- 
tors have been found to account for the variations. The 
spread in metallicity, especially on the low metallicity end, 
for a given i? 2 3 value is due to both to uncertainty/scatter 
in the T e calibration and to additional parameters such as 
the initial mass function (IMF) of the starburst and den- 
sity variations in the nebula (McGaugh 1991). The spread 
in metallicities for a given value of -R23 at the high metal- 
licity end is due primarily to a lack of good T e abundance 
measurements . 

A recent paper by Pilyugin (2000) demonstrates a 
method to remove i?23's density and IMF dependencies 
at the low metallicity end. Pilyugin defines a new pa- 
rameter, p 3 = X 3 - A 23 , where X 3 =log(f([OIII]A4959 + 
A5007)//(H/3)) and X 23 =log( J R 23 ). If the measure of X 23 
is constant for HII regions of similar oxygen abundance, 
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Fig. 3. — This line diagnostic diagram shows the relationship between the [Nil] A6583/Ha line ratio and the [OIII] A5007/H/3 line ratio. 
High metallicity galaxies are found on the lower right where the oxygen line ratio is a good predictor of metallicity. Low mctallicity galaxies 
are found to the upper left where the nitrogen line ratio is a good metallicity indicator. 



then a plot of p 3 vs. X 3 should give a line with a slope of 1 
for objects of the same metallicity. When Pilyugin plotted 
several HII regions in this way, he found the slope is not 
1, implying that X23 can vary for objects with the same 
oxygen abundance. He defined a parameter, Xj*, which is 
equal to the value of X3 when the data in the X3 vs. P3 
plot are extrapolated to p 3 = 0. Doing so he found the 
following relation: 

X; =X 3 - 2.20p 3 - (2) 

For any observed value of X 3 and p 3 one can calculate X 3 . 
Pilyugin goes on to demonstrate a correlation between X 3 
and metallicity which agrees with the T e method to within 
0.1 dex. The p 3 method effectively removes the system- 
atic uncertainties in the R23 method on the low metallicity 
end. As he describes it, the temperature measurement in 
the T e method accounts for the state of the the IMF and 
the geometry factors in the nebula. In this new method, 
P3 replaces the temperature as a descriptor of these influ- 
ences. The correlation Pilyugin found which we will adopt 
for objects on the lower metallicity branch of the R23 re- 
lation is given in by the following equation: 

12 + log(0/H) = 6.35 + I.45X3*. (3) 

Pilyugin (2001a) derives a similar relation for the upper 
metallicity branch. Unfortunately it is of use only in a lim- 
ited metallicity range (~8.1 to ~8.6 dex). In addition, fits 
to the upper branch suffer more from a lack of good cali- 
bration points than from a multi-valued R23 relation. For 
simplicity we will adopt the Edmunds and Pagel (1984) 
fit to the upper branch. This fit, as quoted in Pilyugin 
(2000), can be expressed as follows: 

12 + log(0/H) = 9.57 - 1.381og(i? 23 ). (4) 



With the two methods, R23 and P3, we determine metal- 
licitics for galaxies which possess spectra that contain both 
the strong [OIII] lines AA4959,5007 and the strong [Oil] 
line AA3726,29. When log([NII]A6583/Ha) < -1.3, we use 
p 3 to derive metallicity. When log( [Nil] A6583/Ha) > -1.0 
we use R23 to derive mctallicity. For galaxies with -1.3 < 
log([NII]A6583/Ha) < -1.0 we can not use either method 
as the object is in the turn-around region. 

As an additional consideration, small uncertainties in 
the reddening coefficient, cap, can translate to large er- 
rors in the [OII]A3726 + A3729/H/3 line ratio. Therefore 
metallicities derived with the strong-line method can be 
influenced by weak or noisy H/3 lines. For example, if 
we plot log([OIII]A4959 + A5007/[OII]A3726 + A3729) vs. 
log(i?23) (McGaugh 1991) for our sample of galaxies, we 
find several objects in unphysical locations of the diagram. 
When we plot the same diagram restricting our sample of 
objects to those with an equivalent width of H/3 greater 
than 8 A, the discrepant points are removed. Therefore, 
to ensure good spectral quality and accurate line ratios, 
we have adopted an H/3 equivalent width limit of 8 A for 
the ensuing analysis. 

With the equivalent width limit of H/3 and the nitrogen 
line ratio criteria specified, we calculate metallicities for 59 
additional galaxies, 46 using the R23 method and 13 using 
the P3 method. The metallicities for these galaxies along 
with the oxygen line ratios are given in Tables 2 and 3. 
We also have P3 abundances for the twelve Group I galax- 
ies that have T e metallicities. Seven of the galaxies have 
P3 abundances within 0.1 dex of their T e result. Four are 
within 0.2 dex and one galaxy is highly deviant at 0.4 dex. 
We believe the reason for the deviation is related to the 
evolutionary state of the starburst. This particular galaxy 




Fig. 4. — This line diagnostic diagram relates the [Nil] A6583/Hct line ratio to metallicity. Abundances calculated using the R23 fit arc shown 
as circles, while p3 abundances are shown as triangles. Abundances calculated with the T e method arc shown as diamonds. T e abundance 
data taken from Izotov et al. (1997) are shown as '*'s. We fit a quadratic function to the data (solid line, see text). 



appears to be a highly evolved, low-metallicity starburst. 
Since it is well past the peak of its star formation episode, 
it has an elevated [OII]A3726 + A3729/H/3 line ratio de- 
spite the low overall metal abundance. The Pilyugin fit 
may not track this evolution as his fit is based on the rel- 
atively young starbursts studied by Izotov et al. (1997; 
1998; 1999). Further comparison of the two methods will 
be presented in Paper II. When calibrating the low metal- 
licity end of our line diagnostic diagrams, we will use the 
T e abundances when available and the P3 abundances for 
objects with no T e result. 

3.3. Metallicities from Line-Diagnostic Diagrams 

We use the R23, P3 and T e metallicities presented in 
the previous sections to relate the [OIIIJA5007/H/3 and 
[NII]A6583/Ho! line ratios with metallicity. We choose 
these two line ratios because they are observed in 
nearly all of our ELGs, typically have good signal- 
to-noise ratios, and are fairly insensitive to uncertain- 
ties in reddening corrections. We find that both the 
[OIII]A5007/H/3 and [NII]A6583/Ha line ratios can be 
used as predictors of metallicity. However, they both 
have limitations. A plot of the relationship between 
log([OIII]A5007/H/3) and log([NII]A6583/Ha) is given 
in Figure 3. In this plot, metallicity varies smoothly 
over the distribution of galaxies with low metallicity 
systems in the upper left and high metallicity galax- 
ies in the lower right. At the low metallicity end, the 
[OIIIJA5007/FL3 line remains almost constant for a large 
range of log([NII]A6583/Ho;) values. Specifically, for 



log([NII]A6583/Ha) < -1.2 the [OIII]A5007/H/3 line ra- 
tio is not a good metallicity indicator. Similarly, on the 
high metallicity end the [NII]A6583/Ha line ratio ceases to 
be a good metallicity indicator for log([OIII]A5007/H/3) < 
—0.25 (see also Figure 2). We will use this information to 
decide which line ratio provides the most accurate metal 
abundance for a given galaxy. 

In Figure 4 we plot log([NII]A6583/Ha) vs. metallic- 
ity for the Group I and II data. As the data are sparse 
on the low metallicity end we supplement with T e results 
taken from the literature (Izotov et al. 1997). The diagram 
shows that metallicity increases with nitrogen line strength 
as a smooth single valued function up to the metal rich end 
of the distribution. At log([NII]A6583/Ha) > -0.45 there 
is a sharp upturn in the observed distribution. This corre- 
sponds to the regime where the nitrogen line ratio ceases to 
be a good metallicity indicator, precisely the phenomenon 
seen in Figures 2 and 3. We fit a quadratic function to the 
data points with log([NII]A6583/Ha) < -0.45 and obtain 
the following result: 

12 + log(0/H) = 9.26 + 1.23A + 0.204N 2 , (5) 

N = log ([NII]A6583/Ha). 

The data points have a root-mean-square (RMS) scatter 
about the fit of 0.156 dex. We use Equation 5 to estimate 
the metallicity for galaxies with log([OIII]A5007/H/3) > 
-0.25. 

In Figure 5 we plot log([OIII]A5007/H/3) vs metallicity 
for the Group I and II data. The relation is well defined at 
the high metallicity end. On the low metallicity end, we 
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Fig. 5. — This line diagnostic diagram relates the [OIII] A5007/H/3 line ratio to mctallicity. Abundances calculated using the R23 fit are 
shown as circles, while pz abundances are shown as triangles. Abundances calculated with the T e method are shown as diamonds. T e 
abundance data taken from Izotov et al. (1997) are shown as '*'s. We fit a line to the high metallicity branch (solid line, see text). 



see a scattered clump of galaxies with metallicities ranging 
from 7.5 to 8.2 dex and log([OIII]A5007/H/3) ranging from 
0.3 to 0.8. Again we see that [OIII]A5007/H/3 is not a good 
indicator of metallicity for low metallicity systems as pre- 
dicted by Figure 3. We fit a linear function to all points on 
the upper metallicity branch, where log([NII]A6583/Ha) > 
-1.2. Wc ignore objects in the turn-around region and be- 
low. The result for the fit is as follows: 

12 + log(0/H) = 8.65 - 0.663O x (6) 

O x = log([OIII]A5007/H/3), 

with an RMS of 0.149 dex. We use Equation 6 to estimate 
metallicities for KISS galaxies with log([NII]A6583/Ha) > 
-1.2. 

We use Equations 5 and 6 to compute final metal- 
licity estimates for all KISS ELGs with the necessary 
spectral information available. Included are all objects 
that possess follow-up spectra rated as good or excellent 
quality, that have been classified as starbursting ELGs, 
and for which measurements of both [NII]A6583/Ha and 
[OIII]A5007/H/3 exist. A total of 519 galaxies satisfy 
these criteria. For a number of galaxies, metallicity es- 
timates are computed using a single line ratio, while 
for others both line ratios are used. We combine the 
nitrogen and oxygen metallicity results in the follow- 
ing way. Referring again to Figure 3, for galaxies with 
log([NII]A6583/Ha) < -1.2 we calculate the metallicity 
using only the nitrogen line ratio and Equation 5. For 
objects with log([OIII]A5007/H/3) < -0.25, we use only 
the oxygen line ratio and Equation 6 to calculate metal- 



licity. For galaxies with log([NII]A6583/Ha) > -1.2 and 
log([OIII]A5007/H/3) > -0.25 we calculate metallicities us- 
ing both the nitrogen and oxygen line ratios and Equations 

5 and 6. We take the average of the two results to produce 
a final abundance estimate. When we have estimates of 
the metallicity from both the nitrogen and oxygen line ra- 
tios, we find that the RMS scatter in the difference of the 
two measurements is 0.13, but that the mean difference is 
-0.009 (i.e., consistent with zero difference). 

This process generates metallicity estimates for 519 ho- 
mogeneously observed starbursting galaxies out to a red- 
shift of z=0.095. Since the RMS scatter in Equations 5 and 

6 is 0.156 and 0.149 respectively, we assign a uncertainty of 
0.16 dex to each metallicity measurement. To be conser- 
vative, we use this uncertainty even when averaging two 
metallicity estimates together. Abundance estimates for 
individual galaxies calculated from these empirical meth- 
ods will be tabulated in Gronwall et al. (2002b). Here we 
use the results to identify low metallicity candidates for 
further study. A list of objects with metallicities below 
7.9 is given in Table 4. Several of these objects have been 
observed in detail with the Lick telescope (Paper II) but 
many more warrant high signal-to-noise ratio observations 
yielding T e abundance data. As the objects in Table 4 
are low metallicity, high quality abundance measurements 
of these galaxies can be used to study the primordial he- 
lium abundance (Izotov et al. 1994; 1997) and place con- 
straints on Big Bang nucleosynthesis. In the next section, 
this large sample of data will be used to investigate the 
existence and form of the metallicity-luminosity relation. 
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Metal Abundances of KISS Galaxies I. 



4. THE METALLICITY-LUMINOSITY RELATIONSHIP 

We combine galaxy metallicity estimates with calcula- 
tions of the absolute B magnitudes to investigate the form 
of the metallicity-luminosity relationship. Apparent mag- 
nitudes, corrected for galactic reddening, are measured 
from the imaging portion of the survey data. We adopt 
a Hubble constant of H = 75 km- s _1 -Mpc _1 , and use 
the redshift measured from the follow-up spectrum of each 
galaxy to arrive at the absolute magnitude. 

Our metallicity-luminosity relation is shown in Figure 
6. The general trend is an increase in metallicity with lu- 
minosity over the full magnitude range of the data Mb = 
-21 to -12. Using a bivariatc linear least squares fitting 
technique, we first fit the data with the absolute magni- 
tude as the independent variable, then fit the data with 
the metallicity as the independent variable. Our final lin- 
ear fit is the mean of the two fits, which is described by 
the equation: 

12 + log(0/H) = 3.60 - 0.267Mb. (7) 

(±0.20) (±0.009) 
This is shown as the solid line in Figure 6. The scatter 
about the fit remains constant over the absolute magni- 
tude range of M B = -21 to -16, with an RMS of 0.27 dex. 
The scatter is systematically one sided for Mb > —16, 
implying that there may be a shallower slope at the low 
metallicity end. The small formal errors in the coefficients 
of Equation 7 are somewhat deceiving. In arriving at these 
error bars we assumed an error in the absolute magnitude 
of 0.5 mags, which reflects the formal photometric errors, 
the uncertainty in the Hubble constant, and any possi- 
ble peculiar velocities of the individual galaxies. An error 
in the metallicity of 0.16 dex is assumed, consistent with 
the scatter in the line diagnostic diagrams from which we 
estimate our metallicities. Because we have a large data 
set, these error estimates translate to very small formal 
errors in the slope and intercept of the fits. However, the 
difference in slopes from the direct and inverse fits of the 
bivariate fitting are substantially larger than the quoted 
uncertainty. 

The slope we obtain is affected by two parameters. First, 
the choice of the R23 calibration affects the metallicities 
we derive for our Group II and therefore Group III data. 
Second, the absolute magnitude measurements may suf- 
fer from internal absorption. More massive galaxies are 
likely to have systematically higher extinction than dwarf 
galaxies due to their higher metallicities and higher dust 
content. This would, if left uncorrected, affect the slope of 
the metallicity-luminosity relation. We attempt to address 
both these issues. 

There are several different calibrations relating R23 to 
metallicity at the high metallicity end. We chose the Ed- 
munds and Pagel (1984) fit primarily because of its long 
standing use and simplicity. Pilyugin (2000) offers an up- 
dated i?23 relation which takes into account data not avail- 
able when Edmunds and Pagel made their original result. 
The Pilyugin relation is: 

' 12 + log (O/H) = 9.50 - 1.40 • log{R 23 ), (8) 
which is roughly parallel to but systematically lower than 
the Edmunds and Pagel fit by 0.07 dex. If we use this 
result to calibrate the upper branch of the R23 relation, 
we find the following metallicity-luminosity relation: 

12 + log (O/H) = 3.763 — 0.255Mg. (9) 



The derived slope is slightly shallower than when using 
the Edmunds and Pagel R23 relation. However, it re- 
mains much steeper than previous studies (see below). 
We conclude that while our result is somewhat dependent 
on the -R23 calibration, any reasonable choice for R23 will 
give rise to a significantly steeper slope in the metallicity- 
luminosity relation, compared with the slopes derived for 
dwarf galaxies. 

More problematic is quantifying the absorption caused 
by dust internal to a given galaxy. As massive galaxies 
tend to contain more dust, they suffer from internal ab- 
sorption more than dwarf galaxies. Thus luminosities of 
massive galaxies are likely to be underestimated more than 
luminosities of dwarf galaxies. Correcting the measured 
absolute magnitudes for this internal extinction will make 
the slope of the metallicity-luminosity relation more shal- 
low. As the imaging data do not in general reveal the 
Hubble type or axial ratio of the galaxies, traditional ab- 
sorption corrections are not possible. We attempt a non- 
standard correction in the following way. A plot of B-V 
color vs. Mb for each galaxy shows a trend of increasing 
red colors with increasing luminosity. There is scatter in 
the trend due both to the variations in the stellar popu- 
lations from galaxy to galaxy, as well as to the amount of 
internal reddening and absorption in the galaxy. We fit 
a line to the trend ignoring the most deviant red points 
which are most strongly affected by dust. The fit is taken 
to be the center of the color distribution in the absence of 
internal reddening. We then infer a color excess, E(B-V), 
by measuring the color difference between each galaxy and 
the regression line. We assume that objects redder than 2a 
from the mean regression line, and with cup greater than 
0.25, suffer from internal absorption. The deviant points 
are corrected to the trend line by assuming the following 
reddening law: 

A B = i.0E(B - V) (10) 

While admittedly ad hoc, this absorption correction does 
systematically account for the most heavily extincted 
galaxies in our sample. All of the galaxies that were cor- 
rected have Mb < -16 mags. The result of the correction 
on the metallicity-luminosity relation (using the Edmunds 
and Pagel R23 relation) is plotted in Figure 7. Once again 
a bivariate linear least squares fit is applied and we find the 
following result for the metallicity-luminosity relation: 

12 + log (O/H) = 4.059 - 0.240Mb, (11) 
(±0.17) (±0.006) 

with an RMS of 0.252. 

After the extinction correction is applied to the data, 
we continue to find a slope steeper than those reported in 
the literature. Previous groups have concentrated on the 
low-metallicity, low-luminosity end where T e abundances 
are available. Skillman et al. (1989) used metallicities 
from 20 nearby irregular galaxies, with absolute magni- 
tudes between Mb = -19 to -10.5, to study the metallicity- 
luminosity relation. They calculated metallicities from 
the T e method and distances from Cepheid variables and 
group associations, and found 12 + log(0/H) = 5.50 — 
0.153Mb (dashed line in Figure 6), with an RMS deviation 
of 0.16 dex. Similarly, Richer and McCall (1995) report a 




Fig. 6. — The metallicity-luminosity relationship is plotted in the top panel. The fit to the relationship is given as a solid line while the fit 
derived by Skillman et al. (1989) is shown as a dashed line. The bottom panel shows the residuals to the fit which have an RMS deviation 
of 0.27. 




Fig. 7. — After correcting the luminosities for reddening and adopting the Pagel i?23 relation, we arrive at the metallicity-luminosity 
relationship plotted in the top panel. The fit to the relationship is given as a solid line while the fit derived by Skillman et al. (1989) is shown 
as a dashed line. The bottom panel shows the residuals to the fit which have an RMS deviation of 0.25. 
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metallicity-luminosity relation for 18 nearby dwarf irregu- 
lars of 12 + log(0/H) = (5.67±0.48)-(0.147±0.029)M B . 
Their data range in luminosity from M B = -18 to -10.5 
with a dispersion that increases for M B > -15. Again, 
metallicities are derived from the T e method and distances 
from stellar calibrators. 

Until recently, published fits of the metallicity- 
luminosity relation have not included high luminosity 
galaxies in their samples. However, in looking back to 
data sets on the metallicity of large spirals such as Zarit- 
sky et al. (1994) we find previous evidence for the fact that 
massive galaxies follow a steeper metallicity-luminosity re- 
lation than do dwarf galaxies. This is clear from Fig- 
ure 13 of Zaritsky et al., where the data follow a steeper 
slope than the line included on the plot which has a slope 
similar to Skillman et al. When we combine the Zarit- 
sky et al. data set (39 spirals) with the Skillman et al. 
data, we find a metallicity-luminosity relation given by 
12+log(0/H)= 4.71 - 0.210M B . In another effort to com- 
bine the low luminosity results with high luminosity sys- 
tems, Pilyugin and Ferrini (2000) combine the Richer and 
McCall data set with 13 objects from the Garnett et al. 
(1997) data set and include 17 objects from their own ob- 
servations. These data, ranging in Mb from -21.5 to - 
10.5, have a metallicity-luminosity relation with a slope 
of -0.192. This is evidence that the overall slope may be 
steeper than indicated by the dwarfs alone, though still 
not as steep as our data indicate. One might well ex- 
pect our data to give a steeper slope than the compila- 
tion of Pilyugin and Ferrini. The large spirals used in 
their data set have metallicities derived from several HII 
regions within the disk of each galaxy and averaged to- 
gether. In contrast, our most metal rich galaxies tend 
to be starburst nucleus galaxies, where the emission is 
coming from a large star-forming event in the center of 
the galaxy. Because spirals are known to exhibit radial 
abundance gradients, nuclear starburst galaxies should as 
a group tend to have higher metallicities than those mea- 
sured from disk HII regions, even for galaxies with similar 
luminosities. This is born out in a recent paper by Contini 
et al (2001). They plot a metallicity-luminosity relation 
for a large sample of galaxies. The irregular galaxies in 
their sample follow a metallicity-luminosity relation simi- 
lar to Richer and McCall (1995). The larger HII galaxies 
and UV selected galaxies follow a linear relation with a 
higher slope of -0.173. Displayed on their plot, but not 
included in cither fit, is a large population of starburst 
nucleus galaxies. These galaxies tend to have luminosi- 
ties similar to the UV selected galaxies but tend to have 
higher metallicities. It is clear that the inclusion of the 
starburst nucleus galaxies in a composite fit would lead to 
a significantly higher slope. 

In conclusion, we believe that the steeper slope indicated 
by our data is reasonable. We stress, however, that this re- 
sult is not meant to replace the work of Skillman et al. and 
others, but rather to investigate the metallicity-luminosity 
relationship with a wider sample of galaxy types. In fact 
Figure 7 shows evidence that the overall form of the rela- 
tion may not be a simple linear function, but rather may 
be a higher order polynomial. The slope appears to be- 
come more shallow at the low luminosity end, resembling 
the slopes found by previous studies. Because of the large 



scatter and relative paucity of galaxies at the low lumi- 
nosity end, we do not feel we can adequately justify a 
higher-order fit with the current data set. 

5. DISCUSSION AND CONCLUSIONS 

The metallicity-luminosity relation has significant impli- 
cations for galaxian evolution. It appears to be a contin- 
uous smooth function from high-luminosity massive spi- 
rals to the low-luminosity dwarf galaxies, implying that 
the metallicity-luminosity relation is at work on all mass 
scales and all galaxy types with significant star formation. 
The form of the metallicity-luminosity relation, along with 
its intrinsic scatter, should provide a useful constraint for 
theoretical models of chemical evolution in galaxies. 

We have confirmed the existence of the metallicity- 
luminosity relationship using data from 519 starburst 
galaxies. We found a linear relation between metallicity 
and absolute B magnitude given by: 12 + log(0/H) = 
4.059 — 0.240M B . Metallicities were derived for the large 
sample from secondary metallicity indicators relating the 
strong nebular lines, [OIII]A5007 and [NIIJA6583, to metal- 
licity. We used T e abundances from 12 galaxies, P3 abun- 
dances for 13 galaxies and R23 abundances for 46 galaxies, 
to relate emission-line ratios with metallicity. This study 
uses the largest sample of galaxies to date to construct 
a metallicity-luminosity relation. We find a significantly 
steeper slope than previous results, most likely due to the 
fact that we probe to higher luminosities and include a 
more diverse mix of galaxy types in our sample. The large 
scatter, 0.252 dex, remains roughly constant to an abso- 
lute magnitude of Mb = -16, at which point the scatter 
becomes one sided. This indicates that the slope may be 
more shallow at the low luminosity end, in agreement with 
previous results that have focused on dwarf systems. In- 
terestingly, Richer and McCall have also found the scatter 
in their metallicity-luminosity relation to be one sided for 
objects with Mb > -16. This may be further evidence for 
the need of higher order fits to the relation. We have re- 
sisted the urge to carry out a higher order fit to the data 
due to the large intrinsic scatter in our sample. 

In interpreting the results illustrated in Figure 7, the 
reader should keep in mind the following cautionary note. 
As mentioned above, the abundances we measure for 
the more luminous KISS galaxies are predominantly cen- 
tral/nuclear values, rather than disk HII regions. Since 
spiral galaxies typically exhibit radial abundance gradi- 
ents, our abundances may be biased to higher values com- 
pared to previous studies. Clearly, one needs to con- 
sider the location within the galaxy when quoting abun- 
dances (particularly for large spirals) and investigating the 
metallicity-luminosity relation. While there is no abso- 
lute correct method or convention, we would argue that 
our use of nuclear abundances is no less appropriate than 
using outlying disk HII regions (which will be biased to 
lower abundances). In fact, since dwarf galaxies tend to 
exhibit no substantial abundance gradients, the use of cen- 
tral abundance measurements might well be preferred. 

We have investigated whether this effect could be the 
cause for our steeper metallicity-luminosity relation by 
comparing our results with those from previous studies. 
For example, the large dataset of Zaritsky et al. (1994), 
derived from disk HII regions rather than nuclear star- 



forming regions, exhibits a metallicity-luminosity relation 
very similar to ours at Mb < —18. No bias is evident when 
the two datasets are compared directly. This may be due 
in part to the fact that many spiral galaxies exhibit only 
very shallow abundance gradients (or in some cases none 
at all). We conclude that our use of central abundances is 
not the main cause for our steep metallicity-luminosity re- 
lation, but rather that this is a real phenomenon. We point 
out that evidence for a steeper slope to the metallicity- 
luminosity relation has been available but not pursued 
until Pilyugin and Ferrini (2000) combined the Garnett 
et al. and Richer and McCall data sets to observe the 
overall trend. Their study found a slope steeper than those 
indicated by Skillman et al. and Richer and McCall. Seek- 
ing further confirmation of this result, we combined data 
from Zaritsky et al. (39 spirals) with data from Skillman 
et al. We again found evidence for a steeper slope. The 
metallicity-luminosity relation using this data set is given 
by: 12+log(0/H)= 4.71 - 0.210M B . 

The KISS sample of galaxies includes a diverse mass 
and morphology range, from blue compact dwarfs to giant 
nuclear starburst galaxies. By including the more mas- 
sive galaxies and achieving a large sample, we appear to 
be observing a more general metallicity-luminosity trend 
than is indicated by the dwarf galaxies alone. While the 
metallicity-luminosity relation derived by Skillman et al. 
appears to be a good approximation of the low metallic- 
ity data, extrapolating their relation to higher luminosi- 
ties is probably inappropriate. Efforts to understand the 
physical mechanisms that lead to the observed metallicity- 
luminosity relation will be aided by the constraints set by 
the data presented in this paper. 

The methods used here to calculate coarse metal abun- 



dances allow for a way to quickly identify low metallicity 
candidates for further study. The KISS galaxies believed 
to be low metallicity will be targeted as high priority sys- 
tems for abundance quality spectra in future observing 
runs. Eventually this will yield a large number of low 
metallicity objects that can be used to place constraints 
on the primordial helium abundance of the Universe. The 
coarse abundance methods can also be applied to galaxies 
at higher redshift, where [OIIIJA4363 is difficult to observe. 
Metallicity comparisons between galaxies over a range of 
rcdshifts will help to shed light on the chemical evolution 
of galaxies. 
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Table 1 



Metal abundances from the T e method. 



KISSR a 


log([OII]/H/3) 


logQOIIIJ/H/?) a 12 - 


f- log(0/H) 


49 


0.619 


0.670 


8.00 


85 


0.024 


0.658 


7.50 


396 


0.413 


0.759 


7.90 


666 


-0.310 


0.913 


7.54 


675 


0.145 


0.878 


7.78 


1013 


0.438 


0.829 


7.72 


1194 


0.311 


0.834 


7.95 


1490 


0.179 


0.715 


7.54 


1752 


0.038 


0.678 


7.56 


1778 


0.466 


0.599 


7.93 


1845 


0.393 


0.847 


8.04 


KISSB c 


logaOIIl/H/3) 


log(fOIIIl/H/3) b 12 - 


1- log(0/H) 


23 


0.430 


0.448 


7.50 



a This is the running number for the galaxy as listed in Salzer et al. 2001. 
b [OIII]A4959 + A5007/H/3 

c This is the running number for the galaxy as listed in Salzer et al. 2002a. 



Table 2 



Metal abundances from the R23 method. 



K1SSK, 


ioguunj/up; 


log([Ulllj/H^j a 


12 + log(U/H) 


1 


0.651 


0.249 


8.47 


3 


0.729 


-0.044 


8.47 


68 


0.426 


0.538 


8.48 


91 


0.549 


0.291 


8.55 


109 


0.506 


0.326 


8.57 


110 


0.434 


0.016 


8.78 


117 


0.484 


0.499 


8.48 


124 


0.506 


0.497 


8.46 


145 


0.420 


-0.110 


8.84 


185 


0.355 


-0.093 


8.90 


257 


0.562 


0.433 


8.46 


265 


0.489 


0.389 


8.54 


292 


0.522 


0.605 


8.37 


317 


0.404 


0.637 


8.42 


322 


-0.043 


-0.439 


9.43 


326 


0.631 


0.120 


8.54 


327 


0.580 


0.472 


8.42 


333 


0.430 


-0.188 


8.85 


335 


0.386 


0.358 


8.64 


343 


0.415 


0.824 


8.24 


368 


0.711 


-0.171 


8.52 


398 


0.501 


-0.020 


8.72 


400 


0.371 


-0.239 


8.93 


515 


0.524 


0.067 


8.67 


520 


0.440 


-0.032 


8.79 


577 


0.572 


-0.110 


8.67 


582 


0.759 


0.472 


8.27 


595 


0.005 


-0.803 


9.48 


610 


0.461 


0.202 


8.67 


643 


0.505 


0.441 


8.50 


830 


690 


386 


8 38 


833 


0.566 


0.001 


8.64 


894 


0.614 


-0.044 


8.60 


910 


0.416 


0.063 


8.78 


1016 


0.502 


-0.006 


8.72 


1028 


-0.092 


-0.472 


9.49 


1032 


0.520 


0.295 


8.57 


1055 


0.550 


-0.145 


8.70 


1062 


0.565 


0.345 


8.51 


1402 


0.571 


-0.509 


8.73 


1416 


0.615 


0.325 


8.47 


1424 


0.490 


-0.247 


8.79 


1537 


0.522 


0.510 


8.44 


1885 


0.364 


U.zo 1 


O.io 


1940 


0.626 


0.321 


8.46 


2021 


0.564 


0.084 


8.62 


a [OIII]A4959 + A5007/H/3 










Table 3 






Metal abundances from the p 3 method. 


K1SSR 


log([011j/H/*) 


logQOlllj/H/i) * 


12 + log(0/H) 


96 


0.280 


0.851 


7.91 


97 


0.493 


0.753 


8.05 


105 


0.166 


0.592 


7.65 


120 


0.167 


0.800 


7.80 


272 


0.359 


0.601 


7.85 


310 


-0.097 


0.929 


7.82 


311 


0.248 


0.794 


7.85 


404 


0.459 


0.581 


7.97 


405 


0.154 


0.580 


7.63 


528 


0.507 


0.707 


8.05 


785 


0.287 


0.885 


7.95 


803 


0.600 


0.632 


8.18 


1794 


0.322 


0.736 


7.87 



a [OIII]A4959 + A5007/H/3 



Table 4 



Galaxies with estimated metallicities below 7.9 dex. 



KISSR 


log(IOIIII/H/3) a 


loeflNIII/Ha) 

■"-'oVL if 11K - X J 


12 + lofffO/H) 


p r 

55 


0.808 


T — /■• y a 

-1.654 


t wn 

7.78 


73 


0.623 


-1.531 


7.85 


85 


0.533 


-1.686 


7.76 


105 


0.467 


-1.502 


7.87 


272 


0.476 


-1.528 


7.85 


310 


0.804 


-2.281 


7.51 


311 


0.669 


-1.799 


7.70 


396 


0.634 


-1.515 


7.86 


404 


0.456 


-1.489 


7.88 


405 


0.455 


-1.503 


7.87 


471 


0.676 


-2.340 


7.49 


666 


0.788 


-2.345 


7.49 


675 


0.753 


-1.688 


7.76 


698 


0.844 


-1.577 


7.82 


785 


0.760 


-1.565 


7.83 


799 


0.690 


-1.566 


7.83 


803 


0.507 


-1.505 


7.87 


885 


0.727 


-1.834 


7.68 


986 


0.619 


-2.202 


7.53 


1194 


0.709 


-1.506 


7.87 


1490 


0.590 


-1.606 


7.81 


1752 


0.553 


-1.896 


7.66 


1794 


0.612 


-1.492 


7.87 


KISSB 


log([OIII]/H/3) a 


log(fNIIl/Ha) b 


12 + log(0/H) 


15 


0.607 


-1.556 


7.84 


23 


0.323 


-1.525 


7.85 


35 


0.873 


-1.651 


7.78 


41 


0.786 


-1.620 


7.80 


47 


0.710 


-1.792 


7.71 


53 


0.660 


-1.861 


7.67 


61 


0.656 


-1.586 


7.82 


66 


0.770 


-1.463 


7.89 



a [OIII]A5007/H/3 
b [NII]A6583/Ha 



